Microscopic model for the novel frustrated Cu(II)-spin tetrahedron-based 

Cu4Te50i2X4 (X=C1, Br) systems 
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We present a microscopic study of the electronic and magnetic properties of the newly synthesized 
Cu(II)-spin tetrahedron-based Cu4Te50i2Cl4 compound based on Density Functional calculations 
and on ab initio-deirved effective models. In view of these results, we discuss the origin of the 
observed differences in behavior between this system and the structurally similar and much studied 
Cu2Te205Cl2. Since the Br analog of the title compound has not been synthesized yet, we derive 
the crystal structure of Cu4Te50i2Br4 by geometry optimization in an ab initio molecular dynamics 
calculation and investigate the effect of substituting CI by Br as well as the possible magnetic 
behavior of this system and compare with the recently studied sister compound, Cu2Te205Br2. 

PACS numbers: 71.15.Mb, 75.10.-b, 75.10.Jm 



Introduction.- Frustrated magnetism has gained a lot 
of attention in recent years due to the wealth of new ex- 
otic behavior that arises out of this condition such as spin 
ice and spin liquid phasesi. In the search for new ma- 
terials exhibiting frustrated magnetism, a few years ago 
Johnsson et al.^ synthesized for the first time a family of 
oxohalogenides Cu2Te205X2, X=Br, CI whose structure 
was based on weakly coupled tetrahedra of Cu(II) with 
geometrically frustrated antiferromagnetic (AF) interac- 
tions. These materials have been intensively studied both 
experimentally and theoreticallji^i^i^iSiLSiSiiS. They show 
magnetic ordering with incommensurate wave vectors at 
temperatures T^r = 18 K (CI) and 11 K (Br) and the ob- 
servation of a longitudinal magnon^ in Cu2Te2 05Br2 was 
interpreted as evidence for the proximity of this system 
to a quantum phase transition between antiferromagnet 
and spin liquid behavior. 

The various intratetrahedral and intertetrahedral cou- 
plings and the relative strengths of exchange pathways in 
these compounds have been obtained in detail^ by using 
the electronic structure technique of muffin-tin orbital 
(MTO) based AMTO-downfoldingii. The results pre- 
dicted by this study have been confirmed by subsequent 
neutron diffraction experiments^^ proving the powerful- 
ness of this ab initio Density Functional Theory (DFT) 
based technique in predicting the underlying microscopic 
model of a complex material. 

By changing the subtle ratio between the various inter- 
action paths in these materials, for instance by applying 
pressure or by introducing chemical modifi-cationsiii^ii^, 
one can attempt to drive these systems into quantum 
criticality. Following these ideas, a new oxohalogenide 
Cu4Te50i2Cl4 has been very recently synthesized by 
Takagi et al.^^ which orders antiferromagnetically at Tj^ 
— 13.6 K. This system is structurally similar to the pre- 
viously discussed Cu2Te205Cl2 compound but presents 
some markedly different features. As pointed out by Tak- 



agi et at, the primary structural difference between the 
new Cu4Te50i2Cl4 (which we refer to as Cu-45124(C1) 
following Ref.Ili) and Cu2Te205Cl2 ( Cu-2252(C1) ) is 
the presence of a Te04 complex in the middle of the Cu- 
tetrahedral network in the ab plane (see Fig^. This 
fact led the authors of Ref. 12 to expect an increase in 
the separation between the Cu4 tetrahedra and hence an 
increase in the relative importance of the intratetrahedral 
coupling with respect to the intertetrahedral coupling. 

In the following, we study this proposition within the 
framework of the ATVITO-downfolding technique. In ad- 
dition, motivated by the more anomalous properties of 
Cu2Te205Br2 compared to Cu2Te205Cl2 as reported 
in the literature^i^i^iSiLSiSiifliiiii, we explore the impli- 
cations of substituting CI by Br in the new oxohalo- 
genide Cu-45124(C1). Since the Cu-45124(Br) system has 
not been synthesized yet, we propose the crystal struc- 
ture of Cu4Te50i2Br4 by performing a geometry relax- 
ation in the framework of ab-initio molecular dynamics 
and we analyze its electronic structure by the AMTO- 
downfolding technique. 

Structure.- Both Cu-45124(C1) and Cu-2252(C1) com- 
pounds crystallize in a tetragonal structure. The basic 
structural unit in both systems is the [CuOsCl] distorted 
square (marked in thin lines in Fig. ^) with Cu (marked 
in magenta in Fig. ^ at the center. Groups of four such 
squares share corners, giving rise to [CU4O8CI4] units 
with Cu ions in tetrahedral coordination forming a mag- 
netic cluster of four Cu ions. Cu-45124(C1) has two in- 
equivalent Te atoms, Te(l) and Te(2) (marked in green 
and cyan in Fig.^top panel), while Cu-2252(C1) has only 
one type of Te atoms (marked in cyan in Fig. ^bottom 
panel). The Te(2) atoms in Cu-45124(C1) sit on an equiv- 
alent position to the Te atoms in Cu-2252(C1). Viewing 
the structures along the [001] direction (see Fig.O, they 
show a stacking of Cu4 tetrahedra separated by layers 
of Te(2)-0(3) (Cu-45124(C1)) or Te-0(2) (Cu-2252(C1)) 
units. In the case of Cu-45124(C1), additional Te(l)-0(1) 
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FIG. 1: (Color online) Crystal structure of Cu-45124(C1) (top 
panel) and Cu-2252(C1) (bottom panel) projected on the ab 
plane. In order to emphasize the similarity between the Cu- 
2252(C1) and Cu-45124(C1) structures, we use four unit cells 
of Cu-45124(C1) but show only the section that makes this 
structure analogous to a section of four Cu-2252(C1) unit cells 
showing four connected Cu tetrahedra. Magenta and yellow 
atoms stand for Cu and CI. Two inequivalent Te atoms, Te(l) 
and Te(2) in Cu-45124(C1) are shown in green and cyan colors 
(top panel), while the only one inequivalent Te atom present 
in Cu-2252(C1) is shown in cyan (bottom panel). The brown, 
white and black balls denote 0(1), 0(2) and 0(3) respectively 
for Cu-45124(C1) (top panel) and 0(3), 0(1) and 0(2) for 
Cu-2252(C1)— . Note that every four Cu atoms appearing in 
a square arrangement, due to the projection, actually form 
tetrahedra. We also show the various interaction paths in 
black arrows (see the text for discussion). 



units appear in the same layer as Cu. The relative ori- 
entation of the Cu4 tetrahedra along the [001] direction 
is also different between the two compounds. In the case 
of Cu-2252(C1) the Cu4 tetrahedra show the same ori- 
entation, while for Cu-45124(C1) they alternate between 
successive rows. The latter feature leads to two different 
space group symmetries, P4 for Cu-2252(C1) and PA/n 
for Cu-45124(C1) with elongation of the unit cell in the 
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FIG. 2: (Color online) View of the crystal structure of the Cu- 
45124(01) (top panel) and Cu-2252(C1) (bottom panel) com- 
pounds along the [001] direction. Color scheme is the same 
as in Fig. Q The thick arrows denote the various interaction 
paths. 



ab plane with lattice parameter a = 11.35 A for Cu- 
45124(C1) compared to a = 7.84 A for Cu-2252(C1). The 
unit cell dimensions along the c axis remain compara- 
ble with c = 6.32 A for Cu-2252(C1) and c = 6.33 A 
for Cu-45124(C1). The change in space group defines 
Cu-45124(C1) as centrosymmetric compared to the non- 
centrosymmetric Cu-2252 (CI). 

A crucial difference between the two compounds apart 
from the change in bond lengths, is the relative orienta- 
tion of the Cu-Cl bonds (CI atoms are marked in yellow in 
Fig. ^ among different Cu4 tetrahedra. As has been dis- 
cussed in Ref. |^ the CI atoms play an important role in 
mediating the Cu-Cu interaction in the Cu-2252 systems. 
While for Cu-2252 (CI) the Cu-Cl bonds belonging to dif- 
ferent Cu4 tetrahedra point towards each other (marked 
with red arrows in Fig. ^ bottom panel), in the case of 
Cu-45124(C1), due to the relative shift of the tetrahedra, 
they are oriented parallel to each other (see the red ar- 
rows in Fig. n top panel). This aspect is found to have 
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FIG. 3: LDA band dispersion of Cu-45124(CI) plotted along 
various symmetry directions with F = (0,0,0), X = (7r,0, 0), 
M = (7r,7r,0), Z = (0,0, tt), R = (0,7r,7r) and A = (7r,7r,7r). 
The dominant orbital contributions in various energy ranges 
are shown in small boxes drawn on the right hand side. 



important consequences in the context of hopping inter- 
action pathways as will be discussed later. 

Bandstructure.- Fig. |31 shows the non spin-polarized 
band dispersion of Cu-45124(C1) obtained with the lin- 
ear muffin-tin orbital (LMTO) basisi^ within the frame- 
work of local density approximation (LDA). The bands 
are plotted along the various symmetry directions of the 
tetragonal Brillouin zone. The orbital characters indi- 
cated in the figure are obtained by choosing the local 



coordinate system with the y axis pointing along the Cu- 
0(3) bond and the x axis pointing along the Cu-Cl bond. 
The square planar symmetry of the ligands surrounding 
the Cu^"*" ion sets the Cn-idx2_y2 energy level as the high- 
est energy level. Consistent with the Cu^"*" valency, eight 
bands (there are 8 Cu atoms in the unit cell) dominated 
by Cn-d^^ -y2 character and split off from the rest of the 
bands, span an energy range from « —0.3 eV to 0.2 eV 
with the zero of energy set at the LDA Fermi level. The 
energy bands dominated by other d characters like diy, 
dzxi dyz and d^z^_i span the energy range between « —3 
and — 2 eV overlapping with the 0-p manifold. The Cl-p 
dominated bands appear right above and partly overlap- 
ping the 0-p dominated bands spanning an energy range 
of « 1 eV. These Cl-p dominated bands are separated by 
a gap of « 0.5 eV from the Cu-o?2,2„j,2 dominated bands. 
There is only a negligible contribution of Te(l) and Te(2) 
to the bands crossing the Fermi energy. We note that in 
the low-energy scale, the LDA calculation leads to eight 
almost half-filled bands. Introduction of correlation ef- 
fects within an LDA-I-U treatment are expected to drive 
the system insulating. In what follows though we will 
focus on the ab initio determination of effective one elec- 
tron hopping interactions which are well described within 
LDA and GGA. 

In Fig. ^ we show a comparative study of the various 
partial LDA density of states (DOS) for Cu-45124(C1) 
and Cu-2252(C1). 

While the basic features of the DOS remain the same 
between the two compounds -indicating that the over- 
all nature of the interactions will be similar for both 
systems- there are a few quantitative differences. The 
Cu-d bandwidth at Ep is narrower in Cu-45124(C1) than 
in Cu-2252(C1). The relative proportion of the Cl-p and 
0-p contribution to the bands at E-p is also smaller in 
the case of Cu-45124(C1). The 0-p and Cl-p dominated 
bands, instead of being separated, overlap to a larger 
extent in the case of Cu-45124(C1). Understanding and 
quantifying these differences requires the analysis of the 
bandstructure in terms of a microscopic model. 

Downf aiding and effective model.- A powerful tech- 
nique to construct a low-energy, tight-binding (TB) 
Hamiltonian starting from a complex LDA bandstruc- 
ture is achieved via the ATMTO-downfolding technique. 
It does so by constructing energy dependent, effective 
orbitals by integrating out irrelevant degrees of freedom 
- a method called downfolding. The accuracy of such a 
procedure can be tuned by the choice of N, the number 
of energy points used in the ATVITO calculation. For an 
isolated set of bands, as is the case in the present study, 
these effective orbitals define the Wannier functions cor- 
responding to the Hamiltonian in the downfolded repre- 
sentation. The real space representation of the down- 
folded Hamiltonian in the Wannier function basis gives 
the information of the effective hopping matrix elements. 

For the present study, we construct the massively 
downfolded Hamiltonian by keeping only the Gu-dx2 _y2 
degrees of freedom active and integrating out all the rest. 
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FIG. 4: (Color online) Comparison of the density of states 
between the Cu-45124(C1) (top panel) and Cu-2252(C1) (bot- 
tom panel) compounds. 



The computed, downfolded bands are shown in the top 
panel of Fig. [S] with dotted hnes. With the choice of two 
energy points, Eq and the downfolded bands are in- 
distinguishable from the Cu-dx^ -y2 dominated bands of 
the full LDA calculation shown in solid lines in the top 
panel of Fig. \S\ 

The corresponding Wannier function is plotted in 
Fig- El Two different views of the same orbital are shown. 
The central part has the 3(1^2 _y2 symmetry with the 
choice of the local coordinate system as stated above, 
while the tails are shaped according to Cl-px and 0-px/py 
symmetry demonstrating the hybridization effects. The 
strong pda antibonds are evident in the plot with Cu hy- 
bridization being stronger with CI than with O, a fact 
also evident in the density of states plot, shown in Fig.^] 

The real space representation of the downfolded 
Hamiltonian in the Wannier function basis, Htb = 
— ^j^- tij{c\cj + h.c.) provides the information of the ef- 
fective hopping interaction i^, between the Cu^^ ions 
at sites i and j. The various dominant hopping inter- 
actions are tabulated in Table The notation for the 
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FIG. 5: (Color online) Top panel: Bands obtained with mas- 
sively downfolded Cu d^2_y2 basis (in dotted hnes) compared 
to full LDA band structure (in sohd lines). The iVMTO en- 
ergy points En spanning the region of interest are shown on 
the right-hand side. Bottom panel: The tight-binding bands 
obtained with the hopping interactions shown in Table. |I] (in 
dotted hnes) compared with downfolded bands (in solid lines). 

hoppings are shown in Figs. ^ and |21 While the hoppings 
ti, t2, tx, ta and tr are in-plane hoppings in the plane 
defined by the Cu tetrahedra, tz and tc are out-of-plane 
hoppings. For the sake of consistency, we adopt for Cu- 
45124(C1) the same hopping notation introduced earlier 
for Cu-2252 (CI) in Ref.g 

For comparison, in Table we reproduce the results 
for the Cu-2252(C1) compound from Rcf. H The bond 
lengths corresponding to each hopping clement have been 
also tabulated. For the dominant, intratetrahedral near- 
est neighbor interaction, ti, we observe that while the 
bond length is decreased by only 2.5% for Cu-45124(C1) 
compared to Cu-2252 (CI), ti is reduced by as much as 
22% due to a smaller superexchange path angle in Cu- 
45124(C1) (Z Cu-0(2)-Cu = 105.7° for Cu-45124(C1) and 
Z Cu-0(1)-Cu= 109.8° for Cu-2252(C1)). The intrate- 
trahedral hopping t2 which was weak for Cu-2252(C1) - 
a fact also supported by neutron diffraction^, remains 
weak for Cu-45124(C1). The in-plane intertetrahedral 
hopping tx, remains in magnitude similar to its analog 
in Cu-2252(C1) while other in-plane intertetrahedral hop- 
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FIG. 6: (Color online) Effective orbital corresponding to 
massively downfolded Cu-d^2_y2 calculation viewed in two 
different planes. Plotted are the orbital shapes (constant- 
amplitude surfaces) with lobes of opposite signs colored as 
red and blue. The d^-2_y2 orbital is defined with the choice 
of the local coordinate system with the y axis pointing along 
Cu-0(3) and the x axis pointing along the Cu-Cl bond within 
the square plane. 
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3.147 
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98 
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3.523 
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3.591 







5.539 


12 


4.163 


-10 


ta 


6.180 


15 


6.021 
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td 


7.834 
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8.033 
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8.251 


18 


9.048 


-48 
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5.063 


24 


5.015 


12 


to 


6.332 


-48 


6.320 


-45 



TABLE I: Cu-Cu hopping parameters corresponding to the 
massively downfolded Cu-d^2_j^2 Hamiltonian. The bond 
lengths are in A and the hopping interaction strengths are 
in meV corresponding to hoppings shown in Figs. and |5| 
The values for Cu-2252(C1) have been reproduced from Ref.,2j. 



pings like ta and tr get suppressed. The out-of-plane, in- 
tertetrahedral hopping, remains more or less the same 
as in Cu-2252(C1), while the hopping increases by a 
factor of two. The most remarkable change is observed 
for the diagonal hopping, td, which is reduced to 7 meV 
in Cu-45124(C1) compared to a value of 80 meV in the 
Cu-2252(C1) compound. This reduction however is not 
caused by the elongation of the bond lengths due to the 
insertion of the Te(l)04 group in Cu-45124(C1), as was 
suggested in Ref. [l^ We reveal the origin of this marked 
difference in the following in terms of a detailed analysis 
of the involved hopping paths. 

The tight-binding (TB) bands, constructed out of the 
hopping parameters tabulated in Table ^ are shown on 
the bottom panel of Fig. [S] in comparison to downfolded 
bands. The TB bands compare satisfactorily with the 
downfolded bands. Omission of long-ranged interactions 
such as tc and tr deteriorates the agreement of the TB 
bands with the downfolded bands, proving the essential 
need for inclusion of long-ranged interactions in the cor- 
rect description of this compound. 

Interaction pathways.- It was pointed out in Ref. 
that C\-p degrees of freedom play a crucial role in the 
renormalization process of the effective Cu-Cu hopping. 
Keeping this fact in mind, we carried out downfolding 
calculations where the G\-p degrees of freedom have been 
kept active in addition to Cxi-d^2 _y2 , so as to define a ba- 
sis consisting of Gn-dx2 _y2 and Cl-p. The Cu-Cu hopping 
interactions extracted out of such calculation are tabu- 
lated in Table. ^ For comparison, we show the results 
for Cu-2252(C1) reproduced from Ref. ^. The crucial role 
of hopping paths involving G\-p is evident by comparing 
the hopping interactions between the massively down- 
folded Cu-c?a;2_y2-only calculation and the Cn-dx^^y'^ + 
C\-p calculation. The former includes the renormaliza- 
tion due to Cl-p's while the latter does not. While the 
pattern of renormalization remains essentially the same 
for the intratetrahedral hopping t2 and the intertetra- 
hedral out-of-planc hopping tc, it is quite different for 
intertetrahedral hoppings like tx, ta, td, tr and tz which 
involve pathways via CI atoms belonging to two different 
Cu4 tetrahedra. The most significant change happens for 
the in-plane intertetrahedral diagonal hopping, td- The 
bare hopping strength of td in absence of the renormal- 
ization effect of G\-p is more or less the same between the 
two compounds (8 meV for Cu-2252(C1) and 9 meV for 
Cu-45124(C1)). However, while a large renormalization 
is observed for Cu-2252(C1) when integrating out the Cl- 
p degrees of freedom, such renormalization is practically 
absent in Cu-45124(C1). This difference is caused -as 
pointed out previously- by the different alignment of the 
Cu-Cl bonds belonging to neighboring Cu4 tetrahedra 
which are parallel to each other in Cu-45124(C1) while in 
Cu-2252(C1) they point to each other. This makes the in- 
tertetrahedral Cl-p - Cl-p bonding in Cu-4-5124(Cl) ofppir 
type as opposed to the Cu-2252(Cl) case, where the Cl-p 
- Cl-p bonding was of ppa type. This is nicely demon- 
strated in the Wannier function plot (see Fig.O, where 
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FIG. 7: (Color online) Overlap between Cn-d^2_y2 down- 
folded AHMTOs, placed at two Cu sites situated at in-plane 
intertetrahedral, diagonal positions. Opposite signed lobes of 
the orbitals are colored as blue or magenta and red or orange. 



the effective Cu-dx2_y2 like Wannier orbitals are placed 
at the Cu sites at in-plane diagonal positions. The over- 
lap between the two orbitals provides a rough estimate 
of the strength of the hopping matrix elements. In the 
case of Cu-2252(C1), the C\-p tails from two Cu sites be- 
longing to two different Cu4 tetrahedra overlap to a large 
extend due to direct alignment providing a C\-p-C\-p ppa 
bonding which mediates the Cu-Cu bonding between dif- 
ferent Cu4 tetrahedra. For Cu-45124(C1), in contrast, the 
overlap of the Cl-p tails from different Cu sites belonging 
to two different Cu4 tetrahedra is practically negligible 
due to misalignment of the CI orbitals. 

The nature of the discussed interaction paths plays 
a crucial role in the magnetic properties of this ma- 
terial. Starting from the hopping parameters, i's, the 
exchange integrals, J's, for antiferromagnetic superex- 
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24 
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12 
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TABLE II: TB parameters in meV corresponding to two sets 
of calculations. Set-1: Massively downfolded Cu-d^2_y2, Set- 
2: minimal set consisting of Cl-p and Cu-d^2_y2 degrees of 
freedom (Cu-f CI downfolding) . The numbers for Cu-2252(C1) 
have been reproduced from Ref. 



change paths may be estimated by making use of the 
expression J « AP/U. While this is a valid approach for 
cases like the ti and t2 interaction paths, in general for 
more complicated paths this expression is not anymore 
precise and one has to use more involved estimations of 
the exchange coupling constants. Nevertheless, already 
the knowledge of the hopping parameters gives us the 
clue about the important interaction paths. The drastic 
reduction of the in-plane and to a lesser extent of ta 
and the longer-ranged tr, in Cu-45124(C1) compared to 
Cu-2252(C1), indicate an overall weakening of the interte- 
trahedral coupling in the new compound with respect to 
Cu-2252(C1) and therefore if the system orders at low 
temperatures, the ordering should occur at a lower T^r 
than in Cu-2252(C1), as observed experimentally'''^. The 
spin ordering patterns will be also strongly influenced by 
the change of interaction paths, especially by the near 
absence of the td and reduction of the ta path (we refer 
to the discussion in Ref. [T^ which places the system in 
the limit of weakly coupled tetrahedra. Also the reduc- 
tion by 22% of the ti value implies a smaller intratetra- 
hedron exchange coupling Ji than in Cu-2252(C1). By 
considering J « 4t^/J7, with U=4eV we obtain as ex- 
change coupling constants Ji « 5.8meV = 67K and 
J2 ~ 0.02meV — 0.2K in comparison with the values 
Ji = 2MmeV = 32.9K and J2 = l.58meV = 18AK 
obtained by fitting the susceptibility of a model of in- 
dependent tetrahedra to the experimental data^^. The 
ratio of J2/J1 is largely overestimated in the fitting, pre- 
sumably because of the neglect of the inter-tetrahedral 
interactions. 

Br System. - In an attempt to predict the properties of 
the not yet synthesized Cu4Te50i2Br4 (Cu-45124(Br)) 
and motivated by the discussed proximity to a quan- 
tum critical behavior of the much studied Cu2Te205Br2 
(Cu-2252(Br)), we have investigated the electronic and 
magnetic properties of the ab initio relaxed structure 
Cu-45124(Br) obtained from first principles calculations. 
In order to obtain a theoretical prediction of the Cu- 
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Te(l) 


0.25 


0.25 


0.37660383 


Te(2) 


0.67370804 


0.018727271 


0.87086473 


Cu 


0.75587933 


0.40510264 


0.34817424 


Br 


0.89355065 


0.56680922 


0.32307512 


0(1) 


0.29531785 


0.40609792 


0.23273845 


0(2) 


0.28384154 


0.87238649 


0.36101169 


0(3) 


0.291597 


0.58080843 


0.93893775 



TABLE III: Fractional coordinates obtained by AIMD of the 
relaxed Cu-45124(Br). 



45124(Br) crystal structure, we substituted CI by Br 
in the original Cu-45124(C1) structure and we relaxed 
the volume and internal coordinates performing Car- 
Parrinello ah initio molecular dynamics (AIMD) calcu- 
lationsiS with a projector augmented wave (PAW) basis 
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This procedure has proven to be very suitable 
for predicting reliable new crystal structures^l. For Cu- 
45124(Br) we assume the same tetragonal space group 
P4/n (No. 85) as for Cu-45124(C1). Of the seven atoms 
in the primitive cell, only Te(l) is in WyckofF position 
2c, while all others (Te(2), Cu, Br, 0(1), 0(2), 0(3)) 
are in position Sg. We thus have 19 degrees of freedom, 
but as the AIMD relaxation is done in the conventional 
cell, we need 131 constraints for the 50 atoms in order 
to preserve the symmetry. We verify convergence of our 
structure relaxation not only with the help of the forces 
but we check that each of the 19 independent coordi- 
nates has converged. This is especially important in this 
structure as we find that the relaxation happens in two 
steps: First, the Br atoms introduced into the CI posi- 
tions rearrange, increasing their bond distance to the Cu 
atoms which are the nearest neighbors. This Cu-Br re- 
pulsion makes an adjustment of the 0(1) positions next 
to Br and the 0(2) and 0(3) positions next to Cu neces- 
sary. But as soon as the Br atom has found a relatively 
favorable position, some changes to Cu and 0(l)-0(3) 
coordinates are actually reversed. Thus, the structure 
immediately following the first fast rearrangement would 
have produced quite different interactions strengths than 
the final relaxed structure given in Table Hn As would 
be expected from the different radii of Br and CI atoms, 
we find the largest adjustments in the Br atom positions 
which change by 0.16 A during the relaxation. The other 
changes are 0.05 A for Te(l), 0.04 A for Te(2), 0.04 A 
for Cu, 0.10 A for 0(1), 0.04 A for 0(2) and 0.07 A for 
0(3).^ 

While the volume of the new structure shows a negli- 
gible change with respect to the volume of Cu-45124(C1), 
appreciable changes in bond lengths and angles are ob- 
served. The Cu-0(2) distance which alternates between 
1.94 A and 2.01 A in Cu-45124(C1), becomes 1.91 A and 
2.01 A in Cu-45124(Br). The Cu-0(3) distance is slightly 
smaller at 1.90 A (from 1.91 A in Cu-45124(C1)) . The 
Cu-0(2)-Cu angle changes from 105.7° in Cu-45124(C1) 
to 107.2° in Cu-45124(Br), while the 0-Cu-O angle stays 
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6.332 


-39 
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TABLE IV: Cu-Cu hopping parameters corresponding to the 
massively downfolded Cu-d^2^y2 Hamiltonian. The bond 
lengths are in A and the hopping interaction strengths in meV 
corresponding to hoppings shown in Figs. Q and |21 The num- 
bers for Cu-2252(Br) has been reproduced from Ref. 0. 
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-62 
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19 


68 


11 


34 
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-39 
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-26 



TABLE V: Cu-Cu hopping parameters in meV corresponding 
to two sets of calculations. Set-1: Massively downfolded Cu- 
d^2_y2, Set-2: minimal set consisting of Br-p and Cu-d^2_y2 
degrees of freedom (Cu-fBr downfolding) . The numbers for 
Cu-2252(Br) have been reproduced from Ref. 0. 



nearly constant at 87.0° (87.1° in Cu-45124(C1)). The 
Cu-Cl distance is 2.24 A, the Cu-Br distance 2.42 A. 
Finally, the Cl-O(l) distances alternate between 3.27 A 
and 3.40 A while the Br-O(l) distances are 3.22 A and 
3.55 A. In Table Hill we present the relaxed coordinates 
of Cu-45124(Br). 

We performed AMTO-downfolding for this system and 
in Tables iryl andlVlwe present the bond distance and hop- 
ping values together with those of (Cu-2252(Br)). Both a 
Cn-d^2 _y2 and a C'a-dj.2_y2+JiT-p downfolding were per- 
formed. 

Cu-45124(Br) shows the same trend as Cu-45124(C1) 
regarding the intcrtetrahedral hopping td, namely the 
near absence of Cu-Cu interaction along this path. The 
rest of in-plane intcrtetrahedral hopping paths in Cu- 
45124(Br) are a bit larger than in Cu-45124(C1) but, ex- 
cept for tx, they are smaller than in Cu-2252(Br). From 
the knowledge of the previous systems, a phase transition 
to an ordered state is also to be expected for this system 
at low temperatures. 

An important issue to be mentioned at this point is the 
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value of the intratetrahedron ratio ti in all the com- 
pounds discussed here. Large values of this ratio can be 
related to an enhancement of intratetrahedron frustra- 
tion, what has been already discussed for Cu-2252(Br)^. 
Cu-2252(Br) is found to have a small but nonzero t2 
in comparison to its value for Cu-2252(C1), where the 
t2 hopping path is basically zero, mainly due to the CI 
renormalization. The new set of systems, i.e the synthe- 
sized Cu-45124(C1) and the ah initio computer designed 
Cu-45124(Br) seem to behave in the opposite way. While 
Cu-45124(C1) has a small but nonzero is, Cu-45124(Br) 
relaxes into a structure where the ^2 path is completely 
renormalized to zero by the hybridization with the Br 
ions (see Table in the square planar configuration. 
Though we found an interesting transient structure for 
Cu-45124(Br) with a moderate intratetrahedron is/^i, 
this doesn't seem to be the energetically favored struc- 
ture within the AIMD approach. 

Finally, we would like to note that the value of the 
ah initio calculated hopping parameters is very suscep- 
tible to small changes of distances and angles between 
the atoms. Our AIMD calculations were performed, as 
mentioned previously, within the GGA approximation. 
Consideration of other exchange correlation potentials 
may change slightly the relaxed structure, which could 
be important especially for the intratetrahedron hopping 
paths, where changes of 0.02 to 0.03 A in the distance be- 
tween Cu and 0(2) and of 2.9° in the Cu-0(2)-Cu angle 
are decisive for the variation of the hopping parameters. 

Summary.- To conclude, we have made a comparative 
study between the spin tetrahedron system CusTcsOsCls 
and a newly synthesized compound Cu4Te50i2Cl4 be- 
longing to the same family, in terms of the microscopic 
analysis of the electronic structure. Our study shows that 
although the basic nature of the interactions remains the 
same, there is a drastic reduction of the in-plane interte- 
trahedral diagonal interaction in comparison to the case 
of CusTcsOsCls where this diagonal interaction was esti- 



mated to be nearly as strong as the Cu4 intratetrahedral 
nearest neighbor interaction ti . We show that the origin 
of this reduction is due to subtle changes in the crystal 
structure of Cu4Te50i2Cl4 which causes Cu-Cl bonds be- 
longing to different Cu4 tctrahcdra to align in parallel in 
the new compound rather than pointing towards each 
other as was the case in CusTesOsCls. This reduction of 
the in-plane diagonal hopping in turn increases the im- 
portance of the out-of-plane hopping to the extent that 
some intertetrahedral hoppings along [001] (tc) are even 
about three times stronger than those within the plane. 

In absence of the, yet to be synthesized, Cu4Te50i2Br4 
and motivated by the more anomalous properties ob- 
served in the Br analog to CusTcsOsCls, we have the- 
oretically derived the hypothetical crystal structure of 
Cu4Te50i2Br4 by performing a geometry relaxation in 
the framework of ah initio molecular dynamics. We have 
analyzed the electronic properties of this new system 
within the AMTO downfolding procedure. We observe 
that, while the overall electronic and magnetic behavior 
seems to be similar to its CI sister compound, this com- 
puter designed Br system shows (except for id) a stronger 
in-plane intertetrahedron interaction than the CI system 
-an effect that was also observed in the comparison be- 
tween CusTcsOsBrs and CusTcsOsCls-. We, however 
don't observe any noticeable effect on the intratetrahe- 
dron frustration in the final relaxed Cu4Te50i2Br4 struc- 
ture. 
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